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Abstract 

The period-average rate in the low frequency limit for thermally activated escape 
with periodic driving is derived in a closed analytical form. We define the low 
frequency limit as the one where there is no essential dependence on frequency so 
that the formal limit f2 — » in the appropriate equations can be taken. We develop 
a perturbation theory of the action in the modulation amplitude and obtain a 
cumbersom but closed and tractable formula for arbitrary values of the modulation 
ampitude to noise intensity ratio A/D except a narrow region near the bifurcation 
point and a simple analytical formula for the limiting case of moderately strong 
modulation. The present theory yields analytical description for the retardation of 
the exponential growth of the escape rate enhancement (i.e., transition from a log- 
linear regime to more moderate growth and even reverse behavior). The theory is 
developed for an arbitrary potential with an activation barrier but is exemplified 
by the cases of cubic (metastable) and quartic (bistable) potentials. 
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1 Introduction 



Thermally activated escape over a potential barrier is ubiquitous in physics, 
chemistry and biology (see pQ, [2], [3], [3], [5] and refs. therein). This phe- 
nomenon is important for both quantum p], [2], [6], [7j, [8], [9] and classical 

h, 0, ma, im, na, m, a m mm, ubi, hsi, u eu, m, m 

systems. It can proceed in strong friction (overdamped) , weak friction (un- 
derdamped) and needless to say intermediate regimes (see refs. above). The 
case of thermally activated escape unperturbed by additional external influ- 
ences pioneered by Kramers is exhaustively investigated and by now is a well 
understood phenomenon. 
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However in most physical realizations the thermally activated escape is mod- 
ulated by some external driving. In this case the stationary limit that may 
be a rather good approximation in the absence of external driving is already 
inapplicable. In the presence of the latter the system becomes intrinsically 
non-equilibrium and the problem in known to be notoriously difficult for an- 
alytical treatment [2]. The particular case of periodically driven escape [2], 
[10], [EE], [12], [13], [H], [15], US], [13, [18], [19], [20], [21], [22], [23] is relevant 
among others for chemical physics [26], [21], [25] (where a chemical reaction 
can be influenced by, e.g., laser electric field) and enzymology [2Z] (where an 
enzymatic reaction can be influenced by an oscillating electric field produced 
by the dynamics of protein structure [23] )• According to [17] the periodic driv- 
ing force ""heats up" the system by changing its effective temperature thus 
giving rise to lowering of the activation energy of escape which can be much 
bigger than the real temperature even for comparatively weak fields" . 

Revealing the physical aspects of enzyme action may well become one of the 
most important grounds for application of the periodically modulated ther- 
mally activated escape theory. In support of this point of view it is worthy to 
note that understanding the role of driving at activated escape in biological 
systems is considered by the authors of [T7] as " a fundamentally important and 
most challenging open scientific problem". The reasons for the above state- 
ment follows. The problem of enzyme catalysis is the main unsolved 
interdisciplinary enigma of molecular biophysics, biochemistry and needless 
to say enzymology. Up to now there is no definite and commonly accepted 
understanding of "how does an enzyme work?" (see, e.g., heat controversy 
at a recent conference in the subject issue of Phil. Trans. R. Soc. B (2006) 
361). The idea that dynamical effects may play a crucial role at enzyme ac- 
tion is very popular at present and the concept of the so-called rate promoting 
vibration is a central one in modern enzymology. However the practitioners 
engaged in chemical enzymology traditionally discuss the phenomenon of en- 
zyme catalysis in terms and notions of the transition state theory (as can be 
seen from the materials of the above mentioned conference). The latter is es- 
sentially equilibrium one that is embedded into its cornestone postulate and 
is poorly suited for taking into account dynamical effects in a reaction rate. 
The main tool to study such effects is the Kramers' theory [TJ, [2], [5], [2S]- 
Regretfully at present this theory is much less known and necessitated for 
applications in enzymology than its transition state theory counterpart. 

The thermally activated escape problem at periodic driving in the overdamped 
classical regime was conceptually solved in the papers [UJ, [T2], [13], [H], 
US, US], HB], US], [20], [21], [22], [23]. Two mutually complementary theories 
ffU, [12], [13] and [H], [15], [2D], [21], [22], [23] provide deep insights on the 
behavior of many physical values of interest. They are based on a physical 
idea of the optimal path and are able to provide an imaginable picture of 
the process. The papers [2], [2], [UJ provide description of the escape rate 
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enhancement at weak modulation (the so-called log-linear regime) where the 
change of the activation energy is linear in the modulation amplitude. Most 
important of all, the escape rate enhancement exhibits the replacement of the 
log-linear regime by more moderate growth with the increase of modulation 
amplitude to noise intensity ratio. Such behavior for the intermediate regime 
of moderately strong and moderately fast driving is well described by the 
[TT] , [12] theory that is corroborated by high-precision numerical results. The 
scaling bahavior of the prefactor near the bifurcation point is investigated in 
details [20], [21], [22], [23] ■ In particular the so called adiabatic regime is most 
thorougly investigated [2], [IB], [22], [23] • This regime is defined by the authors 
of HU, [T5], pi], [20], [21], [22], [23] as the limit of slow modulation ft << 1 
where "the driving frequency is small compared to the relaxation rate in the 
absence of fluctuations and the system remains in quasi-equilibrium" and by 
the authors of the [TT], [12], [13] theory as the one that "goes up to driving 
frequencies of the order of the inverse instanton time which is related to the 
curvatures of the potential". In the adiabatic regime the scaling behavior of 
the prefactor for the cubic potential is given by a simple analytical formula 

[22], [23]- 

However the existing literature leaves room for parallel activity for the follow- 
ing reasons. To attract attention of chemists and biochemists to the modulated 
thermally activated escape theory within the Kramers' approach it is neces- 
sary to present its final results in as simple and understandable form as is 
done, e.g., in the transition state theory. On the contrary the results of both 
PH, [12], H3] and [JjE], [IS], [20], [21], [22], [23] theories are presented via in- 
volved notions and values characterizing the optimal action corresponding to 
the minimizing path. Correct making use of these results requires profound 
comprehension of their physical content and mastering in depth the methods 
involved in their deriving. As a matter of fact people engaged in applications 
(whom the author of the present manuscript belongs to) as a rule are con- 
cerned with much more modest objective: how at a given combination from 
the parameter space (noise intensity D, modulation amplitude A, modulation 
frequency Q and characteristics of the static potential U(x)) to evaluate the 
escape rate enhancement in the presence of driving at least for a simple an- 
alytically smooth (i.e., not piecewise) metastable or bistable potential? That 
is why it is desirable to have the period-average escape rate in a closed an- 
alytical form and explicitly expressed only via the parameters D, A, Q and 
characteristics of the static potential U(x) including no other physical values. 
In other words a theory convinient for applications should restrict the physi- 
cal content of the resulting formula only by the notions used at initial setting 
the problem. Besides both [UJ, [12], [13] and [H], [T5], [2JJ, [21], [22], [23] 
theories invoke to rather sophisticated methods such as, e.g., path integrals 
technique. The initial mathematical formulation of the problem is a partial 
differential equation and it seems interesting to see what results can be ob- 
tained among others by means of usual mathematics. These reasons motivate 
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the appearance of the present manuscript. Our aim is to derive by means of 
elementary methods a closed analytical form of the formula for the escape 
rate enhancement in the low frequency limit for arbitrary values of the mod- 
ulation ampitude to noise intensity ratio A/D except a narrow region near 
the bifurcation point. Our approach is not based on a physical idea a priory 
inserted into the theory but rather is a direct purely mathematical treatment 
of the problem. We define the low frequency limit as the one where there is 
no essential dependence on frequency so that the formal limit Q — > in the 
appropriate equations can be taken. At the same time we observe the require- 
ment T K « (where T K is the stationary Kramers' rate) that is necessary 
for the efficient averaging over the period to be possible. The latter means 
that all interesting phenomena related to the so called stochastic resonance 
(taking place at Q = nT K ) are beyond the scope of the present theory. As the 
value IV oc exp [— (U max — U min ) / D] is usually vanishingly small there cer- 
tainly should be a range (perhaps T K « £1 « D) where the contradictory 
requirements IV << Q and Q — > can be reconciled. As this range turns out 
to be more restricted than that of slow modulation Q << 1 we use the term 
low frequency limit instead of adiabatic regime to avoid confusion. 

The manuscript is organized as follows. In Sec. 2 the problem is formulated 
and its solution is argued to be sought by perturbation technique for the ac- 
tion in powers of the modulation amplitude. Sec. 3 and Sec. 4 are devoted to 
the first and second order contributions into action respectively. In Sec. 5 the 
results are combined in a closed form for the escape rate enhancement in the 
low frequency limit. In Sec. 6 the formula is used to obtain plots. In Sec. 7 
a simple analytical formula for the limiting case of moderately strong mod- 
ulation D « A « \J~D is obtained. In Sec. 8 the results are discussed and 
the conclusions are summarized. In the Appendix some technical details are 
presented. 



2 Formulation of the problem 

2.1 Setting the stage 

In this preliminary Sec. we pose the problem and remind some facts on the 
Kramers' theory to introduce designations and notions used further. In the 
Kramers' model a chemical reaction is considered as the escape of a Brownian 
particle from the well of a potential U (x) along the reaction coordinate x with 
x a being the point of the bottom of the well and Xb being the point of its 
top. The problem of interest is to take into account the presence of periodic 
driving with modulation amplitude A and frequency fl. For the driving we 
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adopt without any serious loss of generality the commonly used form 



f(t) = AsiniVLt) 



The results obtained for this simplest case can be directly generalized to any 
arbitrary periodic driving because the latter can be expanded into a Fourier 
series. In the overdamped limit (strong friction case) the Fokker-Planck equa- 
tion (FPE) for the probability distribution function P(x,t) is 



P(x,t) = -F'(x)P(x,t) - 



F(x)+f(t) 



P'(x,t) + DP"{x,t) 



where the dot denotes a derivative in time, the prime denotes a derivative in 
coordinate, F(x) = —U'{x) is the time independent force field and D is the 
so called noise intensity that is actually the ratio of temperature in energetic 
units to the barrier height. In the absence of periodic driving (/(£) = 0) the 
stationary limit of (j2J) is 



= P( x , t) = - F(x)P(x, t) - DP'(x, t) 



(3) 



which can be integrated to yield for the stationary (Kramers') probability 
distribution function Pk(x) the equation 



J = -DP'(x) + F(x)P K (x) 



(4) 



where J is the stationary flux. If we adopt the boundary condition as the 
absorbtion Pk(x c ) = at some point x c (with x c > Xb where Xb is the barrier 
top) then we have 



J f 
Pr(x) = —exp(^—U(x)/D^j I dy exp(l/(y)/ D 



(5) 



For the Kramers' rate (taking into account that N 
have 



K 



J_ 

N 



J 



UJaUJb 

2vr 



-exp 



{U(x b )-U(x a ))/D 



J dx Pr{x) ~ 1) we 



(6) 



where u„ 



U"(x a ) and uib 



U"(x b )\. 
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In the presence of driving the position of the barrier top becomes time depen- 
dent (q b (t) ~ x b — f(t)) and the population of the well is 



N(t) = [ dx P(x,t) (7) 



A convenient operational definition of the reaction rate constant (adopted, 
e.g., in [IT], [12] and [12]) is 

r (i ). 



N(t) 



It should be stressed that the behavior of this value is sensitive to the actual 
choice of the absorbtion boundary x c (t) where we set P(x c (t),t) = 0. The 
requirement we adopt further is that the absorbtion point should be sufficiently 
far from the barrier top 

x c (t) - x b (t) » A (9) 



This requirement is in accordance with that to measure the current well behind 
the boundary argued in 



Let us consider, e.g., the the simplest cubic (metastable) potential U(x) = 
x 2 /2 — x 3 /3 (CP). In this case x a = 0, x b = 1, uj b = 1, co a = 1 and q b {t) w 
1 — Asin(Qt), q a (t) ~ Asin(Qt). By analogy we adopt for the general case 

q b (t) ~x b - Astn(nt); q b (t) = -f(t); q a (t) w x a + Asin(Qt) (10) 



For the quartic (bistable) potential U(x) = —x 2 /2 + x 4 /4 (QP) we also have 
uj b = 1 while x a = —I, x b = and u a = V2- 

Thus taking into account that N(t) ~ 1 and (j2J) we obtain from (jSJ) 

T(t) « -DP'(q b (t),t) + [F(q b (t)) + f(t) + f(t)]P(q b (t),t) (11) 

The rate constant averaged over the period of oscillations T = 2n/Q is 

t+T 

T=- J dsT(s) (12) 

t 
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Our final aim is to calculate the escape rate enhancement 

r 1 T 

A = — ^— J ds [F( qb (s))+f(s)+f(s)}P(q b (s),s)- 



£j J dsP'(q b (s),s) (13) 



K t 

t+T 



To attain this goal we will also need the probability distribution function near 
the bottom of the well that is known to be 1.3(3, E3] 



Here a a (t) is the dispersion that can be identified with the inverse frequency 
of the well a a (t) = l/u a (t). The latter can be evaluated, e.g., for the CP 
u a (t) ~ 1 — f(t). By analogy we adopt for the general case 

u a {t) ^co a - f(t) (15) 



Substituting (TTTT) and (1151) into <HM we obtain 

_ ^ a - Asinjftt) j [Asinjttt) (u a - Asinjyit))] 2 \ 
la ' j ^ V2^D P l 2D j [ib) 



2.2 Form of the action 

We seek the solution of ([2]) in the form 

P(x,t) = P K (x)u{x,t) (17) 

The form ( fTTI) means that 

P(x c ,t) = (18) 



because Pk(x c ) = 0, i.e., adopting this form we have to neglect the possible 
dependence of the absorbtion point x c (t) on time in the presence of driving 
and assume x c (t) = x c where x c is that in the absence of driving. The latter 
may be justified by the requirement ([9]). In the present approach we adopt 
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this approximation without further discussing its validity. Substitution of (JI 
into (|T3|) with taking into account (jSJ) yields for the value of interest 



t+T t+T 

A = - J ds u{q b (s),s) + — J ds exp(-U(q b (s)) / £>) x 



dy exp(u{y)/D) [[f(s) + f (s)}u(q b (s) , s) - Du'(q b (s), s)} (19) 



Qb(s) 



We denote 



(20) 



For the function u(x,t) we obtain the equation 
u(x,t) = Du"(x,t) — 



F(x) + f(t) - 2D^(x) 



u'(x,t) - f(t)V(x)u(x,t) 



(21) 



At A = we must have u(x, t) = 1 for asymptotically large time that suggests 
to seek the solution of (|2"T1) in the form 



u(x,t) = exp[Aa(x,t)} 



(22) 



This form can not be exact because in the limit A — > we obtain P(x, i) 
Pk(x) that can be valid only asymptotically at t — > oo. We denote 



$(x) = 2L>^(x) - 



(23) 



For the function a(x, t) we obtain the equation 
d(ar, t) = Da"(i, t) + DA[a'(x, t)) 2 + 



- Asin(£lt) 



a'(x,t) — ^(x)sin( y Qt) 



(24) 



In the present manuscript we argue the point of view that in the low fre- 
quency limit the function a(x, t) can be sought as a perturbation series in the 
modulation amplitude A 



a(x, t) = ip(x, t) + Ax{x, t) + A 2 n(x, t) + 0(A 3 ) 



(25) 
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The function a(x, t) by its final result in (|T3|) plays the same role as the action 
from the theories [H], [20], [22], [23] and p], [12], p3]. That is why we will 
also use this name. 



The results of the present manuscript testify that at least in the low frequency 
limit a 1 oc 1/D so that u'(x,t) ~ A/D u(x,t). Thus taking into account that 

f(t) oc A and J dy exp(u(y)/ D^j ~ VDexp(^—U(xt,)/D^ (see below) we 
generally have for the escape rate enhancement 



t+T 



A 



T 



ds u(q b (s), s) 



l + OI 



.4 



(26) 



Though at large ratios A/D the term O(^) ~ 1 we will not take it into 
account in the present paper despite of the fact that the method of calcula- 
tion developed below enables us to treat it. This term seems to give minor 
correction and its role will be considered elsewhere. 

Substituting (T25]) in fT22l) . ff2T]) and collecting the terms at powers of A we 
obtain the following system of equations 



<p(x,t) = Dip"(x,t) + $(x)(p'(x,t) - ^(x) sin(Vtt) 



(27) 



X(x, t) = Dx"(x, t) + $(x)x'(x, t) + D \ip'(x, t)] - sin(ttt) (p'(x, t) (28) 



fi(x, t) = D[m"(x, t) + $(x)/i'(x, t) + 2D(p'(x, t) - sin(ttt) x'(x, t) (29) 



etc. From (T5]), ( TlBl) and ( 12T1) we obtain the boundary conditions 



(f'{x c , t) = —sin(Qt); x'(%c t) = 0; 



(30) 



From ( fTBl we obtain 



(p(x a ,t) = sin(Qt); 



0J,, 



X(x a ,t) = -^sin 2 (nt); 



(31) 
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3 First order contribution into action 



3. 1 Equations for the first order contribution into action 



We start from (1271) . Its solution for asymptotically large time can be sought 
in the form 



(p(x,t) = g(x)sin(Qt) + h(x)cos(Qt) 



(32) 



From fl30l) and fl3Tj) we obtain the boundary conditions 



g(x a ) = ; h(x a ) = o 



(33) 
(34) 



Substitution of ( 1321) into ( 1271) yields a system of coupled equations 

Dg"{x) + <$>g'(x) + Vth(x) = V{x) (35) 
Dh"{x) + <kti{x) - Qg(x) = (36) 



This system is rather difficult for analytical treatment. However in the present 
manuscript we restrict ourselves by the low frequency limit Q — > 0. In this case 
the equations fl35l) . fl36|) ) are decoupled 



9 0*0 + p9 fa) = —p- 
h"(x) + ^h'(x) = 



(37) 
(38) 



Taking into account that 
1 



exp 



D 



dr $(r"! 



-exp 



U{s) - U(z) 
D 



(39) 



we obtain the solution of (!37j) . (1381) satisfying (1331) . (|34|) as 



1 1 



exp 



P|-(z)ea;pi 



■TO 



(40) 



10 



h(x) = 0; 



ti(x) = 



We substitute Pr{x) from into (14"U1) and notice that 

1 



/ Gfs exp (^) 



Making use of integration by part and denoting 

V(s)\ 



T (x) = J ds exp 



Eo(x) 



D 



Xc 

E\(x) = J dz ty(z)exp 



dz ^(z)exp 



U(z) 
D 

U(z) 
D 



ds exp 



ds exp 



u(sys 

D ! 



we obtain 



u a D\T Q (x) 



EAx) 



■EAxa) + E 2 (x) 



Making use of (12Uj) and ([SD we obtain 

E±(x) = —exp ( j Tq(x) + J ds exp 



D 



E 2 {x) = exp — r (x) - exp — ] T (x a) 



D 



D 



We recall ( TlTi) and denote 
S(q b (t)) ' 



r (g 6 (i)) 



ds exp 



'U(s) 



F (x a ) 



ds exp 



q b (t) 
D , 



D 



[s - <&(*))- 



(s - x a ) 
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Q n (x) = Ids exp — A(s — x b ) 2 (s — x) 



(51) 



Making use of N2.3.15.1 from [3T] we obtain 

n+l 

Qn(qb(t)) « n! (2A)^— x 
AA 2 sm 2 (ftt)~ 



exp 



(n+1) 



2A Asm(Ot) 



(52) 



where D n (x) is a parabolic cylinder function. Making use of its known prop- 
erties we have 



Qo{q b {t)) 
To (q b (t)) 



2y/X 



erfc —VX Asin(Qt) 
U(x b ) 



V2t 



-exp 



D 



erfc 



ujbAsin(Qt) 
V2D 



Q^q^t)) = ^-=Asin(Vtt)erfc Asin(VLt) 
2v A 

^-exp (-XA 2 sin 2 (ttt)) 
2 A ^ ' 



+ 



(53) 
(54) 



(55) 



where erfc(x) is an additional error function. Making use of (153]) and fl55|) we 
obtain (see Appendix for details) 



(56) 



S(q b (t)) ^ -(q b (t) - x a ) + O(VD) 



and 



qbit) - Xg 
D 



1 



(57) 



The value — is negligibly small compared with Xb D Xa and can be omitted. 
Thus finally we obtain the first order contribution into action as 



<p(Qb(t),t) 



x b — Asin(Tlt) — x a 
D 



sin(Qt) 



(5f 
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4 Second order contribution into action 



4-1 Equations for the second order contribution into action 



From (1281) and the results for the first order contribution into action we obtain 
X(x,t) = Dx"(x,t) + ${x)x'{x,t) + L(x)[l -cos(2ttt)} (59) 

where we denote 

= - 1] (60) 



We seek the solution of ( 159]) in the form 

x (x,t) = v(x) + <fi(x,t) (61) 

where the new functions obey the equations 

Dv"(x) + $(x)i/(x) = -L(x) (62) 
<j)(x, t) = D(f)"(x, t) + <$>(x)<p'(x, t) - L(x)cos{2Qt) (63) 

From (1301) and (l3Tj) we obtain the boundary conditions 

v(x a ) = -||; v'(x c ) = (64) 

<f>(x a , t) = ^cos(2Qt); (f>'(x c , t) = (65) 

4-2 v(x) function 
The solution for the function v (x) is 



v(x) = -^ + Lj dy eXP{ p U 2 ^' ,iyj J dz L(z)PUz)exp(u(z)/D)(QQ) 



{-U{y)/D 



■ '■' a 



p 2 M 



We denote 



/<*) = -/* L(^ P J d s e Xp (Vg) jiy exp (m) (67) 
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After straightforward calculations we obtain 

I(q b (t)) I(x a ) \ 

r (q b (t)) r (x a ) J 

The inner integrals for both I (qb(t)) and I (x a ) mave maximum 
That is why we adopt the following approximations 

I(q b (t))*-L (<&(*)) 

I{x a ) ~ -L(x b )A 

where 

7 , / ?, f j (U(y) 

= y Oz exp I ^— I y as exp I J J dy exp I 

a ?\ / ?\ fU(s)\ f , [U(y)\ 

A = j dz exp I jy- I y as exp I — J J dy exp I — I 

Denoting 

4 r (x a ) r (?&(<)) 

we obtain 

= -§ 

^.5 4>(x,t) function 

We seek the solution of f)63p in the form 

<f>(x,t) = r(x)sin(2Qt) + s(x)cos(2£lt) 

For the new functions we obtain the equations 
Dr"(x) + &r'(x) + 2fis(x) = 
Ds"(x) + $s'(x) - 2fir(x) = L(x) 



•(»(*)) = 4 
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Recalling that we are in the low frequency limit we decouple the equations 
(the limit ft -> in flU}, (j77|)) 

$ 

r"(x) + — r'(x) = (78) 
s"(x) + ia'Or) = ^ (79) 



From ( 1641) . p5j) we obtain the boundary conditions for these equations 

r(x a ) = 0; r\x c ) = (80) 



w2 



s(x a ) = ^; s'(x c ) = (81) 



The solutions satisfying these boundary conditions are 

s{x) = -v(x) (82) 
r(x) = 0; r'(x) = (83) 

Thus we obtain 

/ 
D 



*(©(*)) = ^ (84) 



4-4 x{ x i t) function 



Combining the results we obtain the second order contribution into action 
X(q b (t),t)n-^(l-cos(2nt)) (85) 



To calculate H from (!73|) we need L{x^) and L(qb(t)). Taking into account 
( J60l) the latter means that we need g'(xb) and g' (qb(t)). From (1401) we have 

M-Mri^) (86) 

From (T4T1) and (155]) we obtain 

-Ei(^h) = -Fg(a;6)exp 



—) + 2\ eXP {—) (87) 
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where 



_ , > V2^D (U{x h ) 



2u h 



D 



(88) 



As a result we have 

7T - 2 



7lD 



19) 



Thus 



L(x b ) 



71 



Tl 2 D 



(90) 



From (EH) and fl55l) we obtain 



2a; 



1 



</'(?&(*)) = -ft 1 



D v 
27rAsm(fit) 



7T 



^(g & (t))-^(x 6 ) 



/er/c 2 



uv4sm(£7i) 



'2D 



2uu b VD 



. ( u>bAsin(Qt)\ 1 
er/c [ ^== — ] H 2 ex P 



'2D 



u 2 A 2 sin 2 (tit) 
2D 



■(91) 



Substitution of fl9TT) into fl60l yields L ((fe(£)) that is not written out explicitly 
to save room. 



4-5 A value 



Let us evaluate the integral A given by f[72|) . The maximum of the inner 
integrals takes place at z ~ a^. In this case both inner integrals oc exp f ^Fp^ ) , 
otherwise only one of them is such and the whole integral should be smaller. 
Thus we approximate at z « x b 



As a result we obtain 
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Taking into account that 



To{Xa)ai ^R exp (EM] (95 ) 

u b \ D I 



we obtain 

A ixD 



(96) 



r (x a ) 2u 2 b 

4-6 value 



Now we start rather tedious evaluation of the function B given by (IyTT) . Let 
us first calculate the auxiliary value Go 



n 1, ( U W\ 7a ( U ^\ L 

k3 = J az exp I ^— I J as exp I I y ay exp 



(97) 



The inner integrals have maximum at z ~ x b . In this case we expand the last 
integral into a Taylor series 

f r (u(v)\ / x (U(x b )\ , , \ , N 

y dy exp I — — «(z- x 6 )exp I — x b )exp I — — I (98) 



Then 



e «^ — gi(x fe ) (99) 

2uj b 



Making use of fl55|) we obtain 

(2vrD) 3/2 /f/(x fe )\ 
©o ~ 3 exp — — ^ (100) 



To calculate 9 we decompose it as follows 

e = e + Wi + w 2 + w 3 (ioi) 
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where we denote 



Wi 



Wo 



■>-b 



kbit) 



dy exp (^j^) 



dz exp 



U{z) 



D 



ds exp 



'17(a) 
D 



■>-b 



dz exp 



g b (t) 



■n, 



dy exp 



kbit) 



U(z) 
D 



ds exp 



'U(s) 



D 



dy exp 



'u( y y 



D 



X b 



■''b 



dz exp 



U(z) 
D 



ds exp 



'17(a) 
D 



Qb(t) 



We begin with the itegral in the square brackets from W\ and W3. It 
evaluated as follows 



exp 



/ dy exp (-A(y - x b ) : 



q b (t) 



q b {t) 



J ds exp (— As 2 ) = -^-j=erf (^\^XAsin(Qt) 



-Asin(nt) 



Then W\ can be approximated as 



dz J ds exp 



■U(s) 
D 



x b z 



Altering the order of integration in the double integral we obtain 
'U(s)\ 



Xc Xc 



J dz J ds exp 



■>-b 



ds exp 



'17(a) 
D 



D 



[s - x b 



exp f— Qi(xb) 



Making use of ( 1551) we obtain 

' UbAsin(Qt) 



Wi « 6 erf 



'2D 



At evaluation of the inner integrals in W2 we apply the same arguments as 
those used for deriving (198]) . Then we have 



W2 ~ / exp 

2ojb J 



q b (t) 



'U(z) 
D 



(z - x b ) 



(109) 



After evaluation of the integral we obtain 



W 2 ~ ©o 



eXP [ 2D h 1 



(110) 



For W3 after taking into account (11061) we have the approximation 



W 3 « ^=erf (VXAsin^tfj J dz J ds exp 

q b (t) 



(111) 



For the double integral we have 

7 7 fU(s)\ V^exp(^) ? 

/ dz J ds exp 1^1 « y D 1 I dverfc{V\v) (112) 

q b (t) 2 ^ ' " -Asin(Ut) 



Making use of N 1.5. 1.9 from [32] we obtain 



^ 2 exp (?&L) erf (^M) ( Asin( ^)x 



2ut 



D 



'2D 



erfc 



ui b Asin{VLt) \ \/2txD 



'2D 



TXU) h 



eXP [ 2D h 1 



(113) 



Combining the results we finally obtain 



9 



D 



To (q b (t)) uj 2 



exp 



uj 2 A 2 sin 2 (ntY 
2D 



+ 



7TUJbAsin(Qt) ( ujbAsiniVLt) 
-erf - 



V2nD 



'2D 



(114) 
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5 Result for the escape rate enhancement 



Combining the results and taking into account that T = we finally obtain 
the escape rate enhancement in the low frequency limit 



2tt 



A 



2tt 



exp{ 



A 
15 



x b — Asin(cj)) — x a sin(<f>) — 



2A 2 sin 2 (<\)) 7T-2 1 



D 



4 ' 2nu 2 7r 



/ u 2 A 2 sin 2 U)\ 
eXP (- 2D + 



7iu b Asin((b) ( uj b Asin(d))\ 

. . erf -== — 

V2^D V V2D J 



(U(x b - Asin(cb)) - U(x b )\ 
exp I — ^ ^-^ ^ 1 x 



2oj 

7T 



-exp 



C/ (x;, — Asin((p)) — U(x b ) 
D 



/erfc 2 



oj b Asin(cj>) 
V2D 



x 



2-nAsin((b) , ( uo b AsinU)\ 1 / uJi A 2 sin 2 UY 
erfc [ ^= — ] H 2 ex P — 



2u b VD 
2irAsin((f)) 



'2D 



2D 



2u; b VD 



( u b Asin((py 
erfc 7= | + 



'2D 



—^exp 



u 2 A 2 sin 2 ((f)) 
2D 



2 / u b Asin(4>)\ 

erfc y= — 

V V2D 



(115) 



This formula is the first result of the present manuscript. 



5.1 Plots 



The formula (11151) notwithstanding it looks very cumbersom can be easily 
treated by, e.g., Mathematica. Prior doing it we recall that we want to stay 
within the so called subthreshold driving regime. The latter provides that the 
potential surface always has a minimum and a maximum, i.e., the oscillating 
field is small enough not to distort the physical picture of the chemical reaction 
as the Brownian particle escape from the metastable state. This regime is 
defined by the requirement A < A c where A c = 0.25 for the case of the cubic 
(metastable) potential U(x) = x 2 /2-x 3 /3 (CP) and A c = 2/(3^3) ta 0.4 (see, 
e.g., [18J ) for the case of the quartic (bistable) potential U(x) = —x 2 /2 + 
x 4 /4 (QP). In Fig.l and Fig. 2 the results for the CP and QP respectively 
at relatively large values of the noise intensity are depicted. At the value 
D = 5 ■ 10~ 2 we have A c /D w 5 and A c /D m 8 for the case of the CP and QP 
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respectively. In Fig. 3 and Fig.4 the results for the CP and QP respectively at 
relatively small values of the noise intensity aredepicted. At the value D = 
3 • 10" 3 we have AJD « 83 and AJ D « 130 for the cases of the CP and QP 
respectively. These values limit the horizontal coordinate in the plots. 



5.2 Limit 1 << ^ << -4= 



Let us consider the limiting case of moderately strong modulation D « 
A« \f~D. We denote 



P =1 + ^I 




(116) 



Taking into account the requirement A / yD < < 1 and discarding in ( 11151) the 
terms 0(A/y/D) we obtain 



2tt 



A 



2tt 



expl — 



(xf, — x a ) sin{4>) — Ap sin 2 ((f)) 



;ii7) 



Taking into account A/D » 1 we can evaluate the integral by the steepest 
descent method and obtain a simple formula 



A 



2jnA 



(x b - x a ) /2 - Ap 



expl — 



x b - x a - Ap 



'1181 



where p is the constant given by f 1 1 1 6 j) . For the case of CP we have p ~ 1.632 
while for the QP we have p ~ 2.132. For both of them we have x b — x a = 1. 
In the considered range D « A « \f~D the expression under the square 
root in the denominator can not be zero at physically reasonable values of the 
noise intensity D < 1CT 1 . The formulas fl 1 1 6 [) and (11181) are the second result 
of the present manuscript. 



6 Conclusions 



The results obtained testify that the perturbation expansion for the action a 
in the modulation amplitude A yields a reasonable and convergent expression 
at least in the low frequency limit. In this case we restrict ourselves by the 
second order term O (A/D). In our opinion the corrections from the third and 
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higher order contributions (which are 0(A 2 /D)) will not distort the results 
appreciably except in a narrow region near the bifurcation point A c . In this 
region our results diverge with the scaling laws obtained in J2U], [21], [22], [23]- 
For the regime of weak modulation our results precisely coincide with those 
known from the literature. The formula fll 18j) in the low limit of its validity 
range D << A « \f~D (A small enough for the terms 0(A 2 ) to be neglected) 
with taking into account that — x a = 1 for both the CP and QP yields the 
formula (6.34) from [2] obtained for the case of the QP in the limit of small 
driving frequencies. The value coefficient before the linear in A/D 

term in the exponent of the escape rate enhancement and is actually the so 
called logarithmic susceptibility from the [H] , [15] , [20] , [21] , [22] , [23] theory. 
The value 1 is in agreement with %(0) = lirriQ^o [nQ/sh (7rO)] = 1 obtained 
in [14] for the case of the CP. 

From Fig.l, Fig.2, Fig.3 and Fig. 4 we see that the present theory yields ana- 
lytical description for the retardation of the exponential growth of the escape 
rate enhancement (i.e., transition from a log-linear regime to more moderate 
growth). Moreover Fig.2 exhibits the examples of the reverse behavior when 
the escape rate enhancement attains a maximum at some value of A/D and 
then becomes to decrease with the further increase of A/D. It is worthy to 
note that this phenomenon vividly manifests itself for the QP and is not no- 
ticed for the CP. Regretfully the [II], [12], [13] and [H], p3], [20], [21], [22], 
[23] theories were not exemplified by the case of the QP and our prediction 
can not be directly compared with the results of those theories. 

From the comparison of Fig.l, Fig.2 and Fig.3, Fig.4 respectively we see that at 
a given A/D the periodic driving produces stronger escape rate enhancement 
for the case of CP than that of QP. The latter certainly can not be explained 
by the fact that the barrier height for the CP (1/6) is smaller than that for 
the QP (1/4) because the barrier height does not enter the formulas (j!15p and 
(11181) . This phenomenon can be attributed to the only difference between these 
potentials manifested in the value of the frequency near the bottom of the well 
tu a . For the QP this value (a/2) is higher than that for the CP (1). The latter 
means that the QP goes steeper from the bottom of the well that hinders the 
escape rate enhancement by periodic driving. Thus the precise shape of the 
potential is of utmost importance for the phenomenon of interest. 

We obtain two forms of the resulting formula valid for arbitrary potentials 
with an activation barrier. The formula ( 11151) encompasses the case of arbi- 
trary A/D except a narrow region near the bifurcation point. Its drawback is 
that it is very cumbersome. Nevertheless it can be easily tackled by a com- 
puter with the help of, e.g., Mathematica. Its main merit is that it contains 
only the notions used at initial setting the problem and can be used by peo- 
ple engaged in applications of the theory without reading the rest parts of 
the present manuscript. The formula (11 18ft is valid for the case of moderately 
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strong modulation D « A « \f~D. Its merit is that it has closed and quite 
simple analytical form to be used by practitioners in chemistry and biochem- 
istry for quick by hand estimates. The exponent in this formula explicitly and 
vividly demonstrates how the linear term O (A/D) dominating at weak mod- 
ulation (comparatively small A/D) and responsible for the log-linear regime 
is replaced by that substracted by the term O (A 2 /D) at further increase of 
A/D providing the reardation of the exponential growth of the escape rate 
enhancement. Regretfully we can not directly compare our results with those 
of the [IT] . [12] . [13] theory because the latter is inapplicable for the case of 
slow modulation. 

The formula flllSp is valid in a rather narrow range of the parameters. However 
this range turns out to be relevant for applications in enzymology. The typical 
value of noise intensity at enzymatic reactions is D « 3 • 1CT 3 [27J. In this case 
we have the range of validity for (11181) asl « A/D « 18. The typical values 
of A/D for a particular model suggested in [27] were estimated as A/D m 10 
and find themselves at some stretch within this range. The stringent results 
for this case are given by the formula (j!15p and can be seen in Fig. 3 and Fig.4. 

We conclude that for the particular case of the low frequency limit the aim to 
obtain tractable and convenient formulas describing the escape rate enhance- 
ment by periodic driving seems to be attained. The mathematical tools used 
at their deriving are within the scope of elementary methods. 



7 Appendix 



We denote 

fl =db»7M^)'-db7M^)' (119) 

q b (t) %a 



then 



S(q b (t)) « -{q b {t) - x a ) + R 



(120) 



Our aim here to show that the term R is O ( y Dj then (|56|) will be proven. 
Let us write 



R 



i__d_ 

2Xdx b l 



Inl J ds exp —\{s — Xh) 2 / j ds exp — A(s — Xb) 2 \ (121) 
Kbit) 
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Making use of the substitution rj vA = s — Xb we obtain 

R « — — — /n< J dr exp (~ r2 ) / J dr exp ( — r2 ) J" (122) 



As a result of straightforward calculations we have 

D ( uj 2 h A 2 sin 2 ((j))\ , . 

Wo (%(*)) V 2D J 



Taking into account (1531) we obtain the required result 

R~^f~D (124) 
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Fig. 1. The dependence of the escape rate enhancement on the driving amplitude to 
noise intensity ratio for the case of cubic (metastable) potential U(x) = x 2 /2 — x 3 /3. 
The values of the noise intensity D from the down line to the upper one respectively 
are: 5 • 1(T 2 ; 4 • 1(T 2 ;3 • 1(T 2 ; 2 • 1(T 2 ; 1 • 10" 2 . 
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Fig. 2. The dependence of the escape rate enhancement on the driving amplitude to 
noise intensity ratio for the case of quartic (bistable) potential U (x) = — x 2 /2+x A /4. 
The values of the noise intensity D from the down line to the upper one respectively 



are: 5 • 10" 2 ; 4 • 10" z ;3 • 1(T Z ; 2 • 10 



2. 



-2. 



1 • 10" 
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Fig. 3. The dependence of the escape rate enhancement on the driving amplitude to 
noise intensity ratio for the case of cubic (metastable) potential U(x) = x 2 /2 — x 3 /3. 
The values of the noise intensity D from the down line to the upper one respectively 
are: 3 • 1(T 3 ; 2 • 1(T 3 ;1 • 1(T 3 . 
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Fig. 4. The dependence of the escape rate enhancement on the driving amplitude to 
noise intensity ratio for the case of quartic (bistable) potential U (x) = — x 2 /2+x A /4. 
The values of the noise intensity D from the down line to the upper one respectively 
are: 3 • 1(T 3 ; 2 • 1(T 3 ;1 • 1(T 3 . 



30 



